Field-doping of €50 crystals: A view from theory 
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Abstract. The proposal of using the field-effect for doping organic crystals has raised enormous 
interest. To assess the feasibility of such an approach, we investigate the effect of a strong electric 
field on the electronic structure of Ceo crystals. Calculating the polarization of the molecules and the 
splittings of the molecular levels as a function of the external field, we determine up to what field- 
strengths the electronic structure of Cgo stays essentially unchanged, so that one can speak of field- 
effect doping, in the sense of putting charge carriers into otherwise unchanged states. Beyond these 
field strengths, the electronic structure changes so much, that on can no longer speak of a doped 
system. In addition, we address the question of a metal-insulator transition at integer dopings and 
briefly review proposed mechanisms for explaining an increase of the superconducting transition 
temperature in field-doped Ceo that is intercalated with haloform molecules. 



The doped fuUerenes are materials with very interesting properties. Alkali doped Cgo 
with three alkalis per molecule has, e.g., turned out to be metallic, though close to a 
Mott transition, and superconducting. A problem is, however, that different doping levels 
can only be realized by preparing separate crystals. Moreover, because of the strong 
electronegativity of C^o, no hole-doping has been achieved. In this context the proposal 
of using a field-effect transistor for doping pristine Ceo crystals has raised much interest, 
in particular since such an approach should allow us to continuously change the doping 
by simply changing the voltage applied to the gate electrode. Sadly the reports of such 
field-doping and of spectacular values for the superconducting transition temperatures 
in such devices [1, 2, 3] have been withdrawn [4] after an investigation showed that 
the publications were based on fraudulent data [5]. Nevertheless, it is an open question 
whether field-doping of Ceo crystals could be achieved in principle. In the following we 
address several aspects of this question. 



EFFECT OF A STRONG ELECTRIC FIELD 

Reaching substantial charging (of the order of n electrons per Cgo molecule) in a 
field-effect device requires enormous electric fields. As the induced charge is basically 
restricted to one monolayer of Cgo [6], a rough estimate can be obtained from simple 
electrostatics: For neutrality the charge on the gate must equal that on the monolayer, 
hence the field originating from the gate electrode is given by figate = 2nn/A^Qi, where 
Amoi is the area per Ceo molecule in the monolayer and n is the number of induced 
electron charges per molecule. Thus the external field is about 1 V/A per induced charge, 
corresponding to a voltage drop of about 10 V across the molecule. In such a strong 
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FIGURE 1. Splitting of the molecular levels in the self-consistent multipole field (Z < 2) for a (001) 
monolayer (square lattice) of Cgo molecules oriented such that one of their two-fold axes points in the 
direction of the external electric field (perpendicular to the layer). 



external field the Cgo molecules are strongly polarized. Nevertheless, we find that their 
response is still in the linear regime. Furthermore, in the charged monolayer, the field 
experienced by a molecule is screened by the polarization of the neighboring molecules. 
Taking this into account, we find that the field is reduced by about a factor of two. 
Calculating the splitting of the molecular levels in this screened homogeneous field, we 
find a quadratic Stark effect, with the splitting of both, the tiu and the hu level, becoming 
of the order of the band-width for a field corresponding to an induced charge of three 
to four charges per molecule. This seems to be consistent with the typical doping levels 
that had been reported. 

For a more realistic description of the electrostatics in the field-effect device, we have, 
however, to go beyond considering only a homogeneous field. Given the spherical shape 
of Ceo, the natural approach is via a multipole expansion [7]: We choose an external 
field and the corresponding induced charge per molecule. We determine the multipole 
expansion of the field generated by all other molecules about the molecule centered at 
the origin. Using the linear response of a Ceo molecule to multipole fields (calculated 
ab initio), we determine the new charge distribution on the molecules and repeat the 
procedure until self-consistency is reached. 

Figure 1 shows the splitting of the molecular levels in the self-consistent multipole 
field for different doping levels. While for an external homogeneous field {1=1 multi- 
pole) the splitting is independent of the direction of the field, including the effect of the 
induced charge on the neighboring molecules breaks this symmetry. Surprisingly, the 
asymmetry in the splitting is quite strong, even though the fields that break the symme- 
try are fairly weak compared to the homogeneous field. This is because the multipole 
potentials with even / give rise to a linear Stark effect, which changes sign with the 
external potential and which gives rise to a strong splitting even for weak fields. In addi- 
tion it turns out [7] that the splitting due to the 1 = 1 and 1 = 2 potentials add or subtract, 
depending on the sign of the external field: When inducing electrons they add for the tiu 
level and almost cancel for the hu, while when inducing holes the situation is reversed. 
I.e., when a molecular level is filled, the splitting is substantially enhanced. It reaches 
the order of the band width when inducing about two electrons, or somewhat more than 
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FIGURE 2. Polarization of a Ceo molecule in a homogeneous external electric field of 0.02 a.u. 
(« 1 V/A). The change in charge density compared to the field-free case is indicated by the Ap-isosurface 
at 0.0020/aQ. It turns out that the dipole moment for the /i„ (HOMO) charge density is of the same order 
of magnitude as that of the total charge — but of the opposite sign. 



two holes per molecule. Beyond these fillings the effect of the Stark splitting on the 
electronic structure of the €50 monolayer that carries the induced charge will clearly be 
very large, and one can definitely no longer speak of doping. 



As we have seen, the charge density of the Cgo molecule is strongly polarized in 
an electric field, and one might expect that the main contribution comes from the 
polarization of the highest molecular orbital (HOMO). Calculating the change in the 
HOMO charge density, we find, however, that the dipole moment of the HOMO charge 
density is of the same order of magnitude as that of the total charge - but of the opposite 
sign (see Fig. 2). This surprising result can be understood, e.g., in terms of perturbation 
theory: Expanding the wave function to first order in the external field V = E^z and 
calculating the dipole moment p = ez, find that the leading term is given by a sum 
over the matrix element, squared, with all unperturbed molecular orbitals of different 
parity divided by the characteristic energy denominator. Hence the main contribution 
comes from energetically close-by levels, and the sign of their contribution is determined 
by whether they are energetically above or below the level under consideration. For the 
HOMO in a molecule with a large HOMO-LUMO gap this means that the contributions 
mainly come from the molecular levels below - implying anti screening. We thus see 
that HOMO-antiscreening should be quite general for molecules with large HOMO- 
LUMO gap, and, in fact, it can also be found, e.g., in the series of polyacenes: benzene, 
naphthalene, anthracene, tetracene, and pentacene. 



HOMO-ANTISCREENING 
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FIGURE 3. Gap Eg =E {N- l)-2E{N)+E{N+l) with finite-size correction Eg - Eg{U =0) -U /M 
as calculated by quantum Monte Carlo for the ?i„ and the hu band in the (1 1 l)-plane of the Pa3 structure. 
For the half-fiUed ri„-band the gap opens between f/^ = 0.8 ... 1.2 eV, for the half-fiUed /j«-band between 
f/c = 1.2...1.6eV. 

MOTT TRANSITION 

Since the bands in the fullerenes are narrow, while the Coulomb repulsion between two 
electrons on the same molecule is sizable, the doped fullerenes show effects of strong 
correlation. It is, e.g., only due to orbital degeneracy that A3C60 is metallic and not a 
Mott insulator [8]. In field-doped fullerenes the electrons are restricted to a monolayer 
[6]. Hence the number of nearest neighbors to which an electron can hop is reduced 
and the bands are even more narrow. It is therefore expected that the Mott transition 
occurs at critical values of the Coulomb interaction U below those found in the bulk. To 
determine the transition point, we have performed quantum Monte Carlo calculations [9] 
for a doped (1 1 l)-layer (without Stark splitting) and find that the Mott transition occurs 
between Uc — O.S — 1.2 eV for doping with three electrons, and Uc— 1 .2 — 1 .6 eV for the 
half-filled hu band (see Fig. 3). For integer dopings other than half-filling the transition 
is expected to occur for even smaller values of U [9]. Furthermore, the splitting of 
the molecular levels in the electric field should weaken the effect of the degeneracy 
on the Mott transition and lead to still smaller values of Uc [10]. In particular for the 
hu orbital, the strong electron-phonon coupling might lead to a further reduction of Uc 
[11]. One has, however, to keep in mind that the Coulomb interaction U depends on 
the environment of the molecule. Screening due to the polarization of the neighboring 
molecules is, e.g., responsible for a large reduction of U in the crystal as compared to 
the value for an isolated molecule [12]. Likewise, it is to be expected that for a molecule 
in the monolayer next to the gate dielectric, U might be substantially changed from the 
bulk value. This effect is, however, hard to quantify without knowing the microscopic 
structure of the oxide-Cgo interface. Nevertheless it seems likely that field-doped C^o 
should be insulating at integer fillings. 



ENHANCEMENT OF TRANSITION TEMPERATURE 



In A3C60 the superconducting transition temperature Tc increases with increasing lattice 
constant, i.e., with increasing density of states at the Fermi level [13]. It is therefore 
natural to try the same for field-doped Ceo- The simplest way to increase the distance 
of the molecules in the conducting monolayer is to apply uniaxial stress [14]. An 
alternative approach is the intercalation of the crystal with inert molecules. In fact, for 
field-doped Ceo intercalated with haloform molecules spectacularly increased transition 
temperatures have been reported [3, 4]. A subsequent analysis of the lattice structure of 
these crystals revealed, however, that the lattice is mainly expanded perpendicular to the 
conducting layer, and that the density of states in the doped layer shows no correlation 
with the reported Tc [15]. Therefore, the additional coupling to the vibrations of the 
haloform molecules has been proposed as an alternative explanation of the enhancement 
of the transition temperature [16]. It has, however, turned out that such a coupling is very 
small and, for the two-fold degenerate modes, is even excluded by symmetry [17]. The 
reported enhancement of Tc in haloform intercalated Ceo is therefore not understood. 
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